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ABSTRACT 



This report develops a detailed analysis of a type of 
electrohvdrodynamic power generator which employs an ejector 
and a so-called "fluid flywheel" as essential components. 

The medium is steam containing electrically charged water 
droplets . 

The analysis takes into account the experimentally estab- 
lished fact that the maximum strength of the electrical field 
that can be sustained at incipient breakdown at the most criti- 
cal location is proportional to the fluid density at that loca- 
tion. It is shown that as a consequence of this fact, the 
electrical output can be maximized by designing the primary 
jet for an exit Mach number of 0.71. 

Estimates are made of the pump work required, of mixing 
losses in the ejector and of friction and secondary flow 
losses. The mathematical analysis is reduced to a fully non- 
dimensional form and the key dimensionless parameters that 
govern performance are clearly identified. A preliminary 
estimate is made of the numerical values of these parameters 
and the overall performance of the system is estimated on this 
basis . 

Unfortunately, the results so obtained are very pessimistic. 
They indicate that even at 100 atmospheres pressure, electrical 
breakdown so severely limits power output, that it is probably 
insufficient to cover the demand for pump power and for power 
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to overcome the various losses involved, let alone provide any 
useful net output. The only hopeful note is that the analysis 
so clearly pinpoints the problem that it might in the end also 
suggest the means for surmounting it. 

A modified design where the electrical conversion section 
is placed after the ejector is analysed in the Appendix. This 
modification yields less pessimistic results, but still leaves 
doubt as to the viability of this type of power generator. 
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1. Introduction 



The purpose of this analytical report is to develop a mathe- 
matical model of a proposed type of electrohydrodynamic (EHD) power 
generator which embodies an ejector as an essential component. The 
overall configuration is approximately toroidal in form in order 
to help conserve the angular momentum of the working fluid which 
circulates through it. For this reason the inventor* of this con- 
figuration has dubbed it a "fluid flywheel". 

Fig. 1.1 is a schematic diagram of the device. Superheated 
steam at stagnation pressure p Q and stagnation temperature T q 
is supplied to a nozzle which discharges into the ejector mixing 
region at station 1. A secondary fluid enters the mixing region at 
station 2. The resulting mixture leaves the ejector at station 3. 
The ejector action creates a pressure rise across the mixing sec- 
tion which is just sufficient to offset the effect of the power out- 
put and of the various flow losses that occur in the rest of the 
circuit . 

The mixture of liquid and vapor which leaves station 3 is cooled 
during its return circuit sufficiently to condense and separate out 
a mass flow of liquid equal to the mass flow of steam which passes 
through the primary jet. 

The liquid withdrawn from the condenser is pumped back into the 
boiler and superheater where it is heated back to the initial stag- 
nation state 0. The primary fluid thus undergoes a Rankine cycle 
where the turbine and condenser have been replaced by an EHD 
generator. 

* See note in Section 9, References. 
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Fig. 1 . 1 SCHEMATIC OF POWER GENERATOR 
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Conditions in the primary nozzle are controlled in such a 
way as to cause fine droplets of moisture to condense here and 
electrodes are provided which electrically charge the liquid drop- 
lets. These electric charges are removed from the flow by a grid 
type collector located at the ejector exit, station 3. The moving 
charges constitute an electric current which sets up a correspond- 
ing electric field. The electric forces acting on the charged 
particles tend to oppose the motion. However, the size and charge 
of the particles is controlled in such a way that they have negli- 
gible mobility with respect to the surrounding fluid. Thus the 
moving gas does work on the fluid particles in moving them down- 
stream against the resistance of the electrical forces. In per- 
forming this work, the gas stream undergoes a corresponding de- 
crease in enthalpy. The work done by the gas against the elec- 
trical forces sets up an electrical potential difference between 
the charging electrode and the collector grid. If these two 
terminals be connected by an external electrical circuit, that 
circuit will therefore deliver useful electrical power. The 
only turbo device required is a small pump. 

The question that arises is whether such a hypothetical scheme 
can be made practically effective. There are various formidable 
obstacles to overcome, including the achievement of negligible 
mobility. It is therefore important to have a complete quan- 
titative model of all cycle components to help guide research 
and development efforts. 
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Such a model has now been developed for the fluid flywheel 
concept and is presented in this report. The model shows that the 
performance that can be attained is governed by just a few basic 
dimensionless parameters. It identifies these parameters and 
provides some basis for estimating their probable magnitudes. 

By identifying the key parameters and showing their effect on 
overall performance, this analysis permits current experimental 
research efforts to be redirected along more fruitful lines. 

This report is the result of research sponsored by the Depart- 
ment of Energy. 
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2. Basic Electrical Performance 



It is convenient in the present analysis to introduce a 
parameter, a which is used to distinguish between the two dis- 
tinct cases of positively and negatively charged particles. 
Specifically, we set a = +1 for the case of positively charged 
particles, and we set a - -1 for the case of negatively charged 
particles. It is then appropriate to denote the electrical charge 
per unit mass by the product aq where q has the units of 
coulombs/kg and is always positive by definition. 

The electrical working section of the generator coincides 
with the mixing region of the ejector and, assuming polar sym- 
metry, the distribution of q is essentially two dimensional 
over this region. However, our purpose in the present section is 
to develop certain key dimensionless parameters which depend 
primarily on how q varies with respect to the axial coordinate 
z. These parameters are not sensitive to variation with respect 
to the radial coordinate r. Under these circums tances it is 
permissible to simplify the analysis in the following way: At 

each axial station z that lies within the mixing region, we 
subdivide the total cross-sectional area A into an inner circular 
region of area A' and an outer concentric annulus of area A" 
such that 

A' + A" = A = A^ = constant along channel (2.1) 

While fluid density P (kg/m) is treated as essentially uniform 
over the entire mixing region, the axial velocity (m/sec) is 
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assumed to have the value V' over the inner area A' and the 



value V" = V ? over the outer annular area A". Similarly, charge 
density (coulomb/kg) is assumed to have the value q' over the 
inner area A' and the value q" = 0 over the outer annular 
area A". Moreover, the quantities A', V' and q' are regarded 
as functions of coordinate z which vary monotonically from their 
initial values A^, v i ' at z = 0 to their final values 

A 3 , V 3 , q^ at z = 1 . The exact form of these functions is not 
of immediate importance in the present development. 

It is useful to introduce the following auxiliary definitions 
at this point on 

l = length of working section, m 

V = magnitude of voltage change across working section, 
positive by definition, volts 

e = = permittivity of free space 



g = — = dimensionless charge/mass ratio 



Poisson's equation and its first two integrals can now be 
represented in the following one dimensional version: 



8.854 x 10 12 farad/m 



<t> = electric potential, volts 



f = = dimensionless electric potential 



2 

? = j = dimensionless axial coordinate 



q' 




( 2 . 2 ) 



2-2 



(2.3) 







0 0 



(2.4) 



In these equations, the quantity ^ is the negative of the 

dz 

electrical field strength (volts/m) . This quantity attains its 
maximum magnitude at station z = Q, which is the critical loca- 
tion at which electrical breakdown first occurs. 

It has been established by experiment, Ref. (1) , that over 
a certain range of pressures, the breakdown strength of a given 
medium is well approximated by the following expression 




+ CgRo-j. 



(2.5) 



where R ( joules/kg°K) is the gas constant of the medium and where 

2 

the breakdown constants C (volts/m) and C„ (m 0 K/coulomb) are 

o a 

characteris tic properties of the medium which can be measured ex- 
perimentally . 

The evaluation of Eqs (2.3) and (2.4) at the exit station 
5=1 is of particular importance. In this connection it is conven- 
ient to introduce the following notation: 



1 



= / g(Od5 



( 2 . 6 ) 



0 
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g(?)d?dC 



(2.7) 



ip ~ — 



0 0 



Notice that and ip^ are dimensionless coefficients that 

depend only on the form of the dimensionless charge/mass function 

g(0 . 

The length & of the working section of an ideal EHD generator 
having a given current flow i (amps) should be so chosen as to 
maximum the potential difference V (volts) and hence the electri- 
cal power P g (watts) which is generated. This means that the 
dimensionless potential f should attain its maximum value at 
the exit station, or that 



= 0 at £ = 1 



It also follows that 




( 2 . 8 ) 



( f - f 3 ) = a at e = 1 



(2.9) 



Upon substituting Eqs (2.5) through (2.9) into Eqs (2.3) and 
(2.4) we can readily extract a pair of useful relations from the 
result. In this connection it is also convenient to introduce 
the auxiliary notation 



1 
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Note that a is dimensionless. This quantity is normally very small 
compared with unity. 

Using this notation and procedure, we may now summarize Eqs 
(2.3) and (2.4) in the form 



*1 



l = 



(1+oQ 

*1 



£C b r 



( 2 . 11 ) 



V 




(l+cOCgRp^ 



( 2 . 12 ) 



The current flow i (amps) through the channel may be written 

i = aq 1 P 1 V 1 A 1 (2.13) 

The gross electrical power output P e (watts) is now given 
by the simple expression 

P e = [oV]i (2.14) 

In order to reduce our final result into a more useful form 
we also write 

V 1 = a l M l (2.15) 

where 

= axial velocity to nozzle exit, m/sec 
a^ = sonic velocity at nozzle exit, m/sec 
= Mach number at nozzle exit 
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Upon combining the last five relations and rearranging, we 



obtain the important result 

P e " ' ^)< 1+ “> 2eC l R2 » l 2a 1 M i A l <2 ' 16) 

For the purpose of non-dimensionalizing this result, it is 
useful to introduce the following additional dimensional reference 
parameters 

a Q = sonic velocity corresponding to stagnation conditions 
at inlet to primary jet, m/sec 

P = density corresponding to the above stagnation conditions, 
kg/m 3 

A = throat area of primary jet, m 

2 2 2 

Now dividing Eq (2.16) through by eCgR P 0 ^ Q A t 
gives 



? 2 2 

£C B R Po a o A t; 





(2.17) 



Let us denote the bracketed quantity on the right in the 
following way 

| (1 + “ ,2 fe) 2 ft) Mi &)| = (U “ )2 *3 ''"l 1 <2 - 181 

where ii> 3 is a certain normalizing constant which will be defined 
presently . 
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In order to analyze the function defined by Eq (2.18), we treat 
the expansion through the primary nozzle in the usual way as an 
isentropic process of constant total enthalpy. The medium is 
treated as a perfect gas for which y , the ratio of specific heats, 
is a constant. 

The relations which govern such a nozzle flow are well known 
and will therefore be used freely here without detailed deriva- 
tion. Note that there are two distinct expressions for the func- 
tion (M^) depending on whether is subsonic or supersonic. 

For the subsonic case, the throat of the nozzle coincides with the 
exit station so that (A^/A^) equals unity. In the supersonic 
case (A^/A^) becomes a function of . 

When these details are properly worked out, the following 
solution is obtained. 



For M x < 1 
<J> 3 F (M 1 ) 



M, 




(Y+3) 

2 (y-1) 



(2.19) 



For > 1 

i|» 3 F(M 1 ) » 



(Y + l) 
-2 ( Y~l) 



2 

Y + l 



1 



+ 





1 

(Y-1) 



( 2 . 20 ) 
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By differentiating Eq (2.19), it can be shown that the function 



passes through its maximum value precisely at 



M. = — =0.71 

iopc ^ 



( 2 . 21 ) 



Note, incidentally, that this result happens to be independent 
of y : 

2 

On the other hand the function (1 + a) does not 

necessarily pass through its maximum at precisely this same Mach 

number because, as may be seen from Eq (2.10), the quantity a 

is itself a function of density p at station 1, and hence of 

. Nevertheless, a is very small compared with unity, and its 

variation will therefore have an effect on M. , S o small that 

lopt 

it can safely be neglected. 

As a normalizing condition we stipulate that the maximum 
value of the function F(M^) shall be unity. Thus 



F = F (M. J = 1 
max lopt 



( 2 , 22 ) 



This constraint now fixes the magnitude of the normalizing con- 
stant ip^ which works out to be 



V 



Y + 3 



(Y + 3) 
2(y-1) 



(2.23) 
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For example, for y = 1.3 this formula gives 



iJj 3 = 0.4211 



(2.24) 



The dimensionless coefficients \p 2 and ^ and the factor 

2 

(1 + ct) can next be combined into an overall dimensionless coef- 
ficient ^4 according to the relation 






(1+cO 2 



(2.25) 



With this notation, the electrical power output can be expres- 
sed in the following form, namely. 



\ eC B 2R2c 0 



2— = * 4 F( Ml ) 
o A t/ 



(2.26) 



where the function F(M^) is defined by Eqs (2.19), (2.20) 

and (2.23) . 

The function F(M^) is listed in Table 2.1 and plotted in 
Fig 2.1. These results are for y = 1.3 which is a typical value 
often used for steam. Notice that F(M^) reaches its peak value 
of unity at an optimum exit Mach number of approximately 0.71. 
Observe that appreciable deviations from this optimum exit condi- 
tion cause heavy performance losses. 

Under these conditions, the subsequent analysis will be re- 
stricted specifically to the optimum exit Mach number M^ =0.71 
as derived above. Consequently, Eq (2.25) can then be simplified 
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TABLE 2 . 1 



Y = 1.3 



0.0 


F ( M 1 ) 


0.0 


0.20 




0.4550 


0.40 




0.8014 


0.60 




0.9774 


0.71 




1.0000 


0 . 80 




0.9849 


1.00 




0.8722 


1.20 




0 . 7241 


1.40 




0.5886 


1.60 




0.4704 


1.80 




0.3711 


2.00 




0.2901 


2.20 




0.2253 


2.40 




0 . 1744 


2.60 




0.1347 


2.80 




0.1041 


3.00 




0.0805 


3.20 




0.0625 


3.40 




0.0486 


3.60 




0.0330 


3.80 




0.0298 


4.00 




0.0235 


4.20 




0 .0186 


4.40 




0.0148 


4.60 




0.0119 


4.80 




0.0095 


5.00 




0.0077 
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M, ► 



Fig. 2.1 F(M^) versus 
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to the form 






e 




a o A t/ 



^4 



(2.27) 



In order to get the greatest benefit from Eq (2.27) we should 
know the value of the power coefficient ^ . This coefficient 
is defined by Eqs (2.6) , (2.7) , (2.23) and (2.25) . It may be seen 

from Eqs (2.6) and (2.7) that in order to calculate the quantities 
and 1^2 that ultimately fix the initially unknown coefficient 

we must know the exact dimensionless charge/mass distribution 
as fixed by the function g(?) . Unfortunately this function re- 
mains unknown so that an exact solution is not possible within the 
scope of the present simplified analysis. 

Nevertheless, we do have some approximate information about 
this function which is adequate to provide a first order estimate 
of Let us call this . In order to establish such an 

estimate or reference value it is only necessary to assume a 
qualitatively correct and analytically convenient form for this 
function . 

Note that 

g = 1 at 5 = 0 (2.28) 

9 9 3 q ] _ at Q = 1 (2.29) 



Moreover g(£) varies monotonically between these two end 
values and it should also satisfy the condition, at least 
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approximately, that 




(2.30) 



It is also timely to note here that the current i evaluated 
at the two stations 1 and 3 gives the relation 



q l P l A l V l q 3 p 3 A 3 V 3 



(2.31) 



From this we readily infer that 



= ^3 _ / p l A l V l\ 
3 ' \ P 3 A 3 V 3/ 



X 



(2.32) 



where x is the dimensionless mass flow ratio of the ejector. 

A simple analytical form for the reference function g(S) 
that satisfies the foregoing constraints is 



gU) = x + (1- x)(l-<;) 2 (2 

Upon substituting this into Eqs (2.6) and (2.7) and inte- 
grating we obtain the expressions 



'P 



1 



1 

3 



(l+2x) 



(2.34) 



P 2 = J ( 1+x > 



(2.35) 
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Then using Eq (2.25) we obtain the desired reference value as 



^40 



| i|> 3 (1+a) 2 



(l+5x) 
(l+2x) 2 



(2.36) 



where constants ip^ and a are still defined by Eqs (2.23) and 
(2.10) as before. 

Of course is only a reference, not the true value ip^ 

itself. However, we may relate these two quantities as follows. 



^4 = ^40 




(1+a) 



2 



( 1+5 x ) 

( l+2x) 2 



(2.37) 



Here y is regarded as an empirical correction factor of order 
unity. It is to be expected that whereas ip^ itself may tend to 
vary over a certain range, y should be more nearly constant. 

We may now combine Eqs (2.27) and (2.37) in the form 



\ 


| <IMi+c0 2 


( 2„2 2 , u 


R p a A, / 
\ B K o o t 





( 1+5 x) 
(l+2x) 2 



(2.38) 



The quantities ip^, a and x which appear in Eq (2.38) are 
defined by Eqs (2.23), (2.10) and (2.32), respectively. The coef- 

ficient y must be determined by experiment or by more advanced 
analysis. However, for purposes of order of magnitude estimates 
it is permissible to take 



y = 1 



(2.39) 
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3. Gross Electrical Power: Change of Dimensional Base 



In analyzing the mixing and friction losses in the ejector 
and the friction and secondary losses in other parts of the sys- 
tem, it will prove useful to employ as dimensional reference 
quantities the fluid density P = P]_ = P2 = P3 the e j ector ' 

the exit area from the ejector and the average axial 

velocity at this station. Thus the reference mass flow rate 

becomes (kg/sec), the reference kinetic energy per unit 

2 

mass becomes V^/2 (joule/kg) and the reference power becomes 
(p^A^V^) (V^/2) (watts). In this connection it is appropriate to 

define a dimensionless gross electrical power coefficient c q 
in the following way, namely. 



It is also necessary to establish suitable links between the 
above reference quantities and certain corresponding quantities 
at station 1, the exit from the primary jet. In this connection 
recall that 



P 



e 



(3.1) 



c 




e 




flow ratio 



(3.2) 



We also elect to introduce the following definition 
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(3.3) 



V l\ 

~J = Y = dimensionless velocity ratio 



As before it will be convenient to write 



V 1 = a l M l 



(3.4) 



Moreover we can substitute for P g in Eq (3.1) its expression 
as previously developed in Eq (2.38). 

When the foregoing relations are combined and rearranged, the 
following result is obtained 




(3.5) 



This suggest the utility of defining the following auxiliary 
quantities, namely, 




Recall that the quantities involved in Eq (3.6) all correspond 
to the value M ]_ opt = 1/^2. Hence ^ turns out to be a constant 
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which depends only on y and a. When the details are worked out 
we find that 



Since the medium is treated as a perfect gas, the crucially 
important parameter defined by Eq (3.7) can be further reduced to 
the form 



This result shows the importance of maintaining inlet stagna- 
tion pressure p Q as high as possible and inlet stagnation tempera- 
ture T q as low as possible. However, if condensation is used to 
produce droplets, T q is restricted to that range of temperatures 
which produces optimum droplet size. One way around this restriction 
is as follows: Let the carrier fluid be changed from steam to air 

(or to some other gas which is noncondensible at ordinary tempera- 
tures) . Let charged liquid droplets or solid particles of proper 
size be injected into the carrier fluid. This permits T q to be 
sharply reduced whereupon 3 is correspondingly increased. Of 
course this alternative may involve its own characteristic draw- 
backs as well, but these will not be discussed further at this 
point. 




(3.8) 



For example, for y =1.3 this gives 



4 > = 1.32(1 + a) 2 



(3.9) 




(3.10) 
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While the introduction of the separate quantities y, and 8 
follows naturally from Eq (3.5), we note that these quantities 
always occur in the form of the product y^g . Hence it will 
simplify the subsequent development to set 



u = y^8 

whereupon Eq (3.5) reduces to 



(3.11) 



,, (1 + 5 X > 2 

c y *- 2 xy 

( 1 + 2 X ) 



(3.12) 



In the subsequent analysis of ejector performance, it will be 
shown how the two fundamental parameters x and y govern pressure 
rise through the ejector, mixing losses and so on. Eq (3.12) 
therefore ties in gross electrical power output with these other 
effects on the same consistent overall basis. In particular, it 
is important that the gross electrical power coefficient c e as 



expressed by this equation should exceed the sum of all the 
losses in the system by as wide a margin as possible. 
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4. 



Ejector Performance 



Recall the notation introduced earlier, namely, 




x 



mass flow ratio 



(4.1) 



= y = velocity ratio 

As a basis of comparison with the actual physical ejector, 
consider a hypothetical ideal ejector which satisfies the following 
constraints. Firstly, the fluid density is uniform throughout 
the ejector so that p^ = p 2 = P 3 = p = constant. Secondly, pres- 
sure is uniform across the entire inlet so that p 2 = p . Thirdly, 
the velocities are purely axial and are uniform across sections 
1, 2 and 3. Fourthly, wall friction is negligible. 

Imagine this ideal ejector to be in operation with all elec- 
trical circuits turned off. The ideal pressure (P 3 * ~ P-j_) that 
would occur under these conditions may be expressed by means of 
the dimensionless coefficient 




(p. 



- p x > 



I pV 3 2 




ideal pressure rise coefficient 
(for given values of x and y) 



(4.3) 



Of course the real physical unit operating under realistic 
conditions with electrical power output connected to a load, will 
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produce some smaller pressure rise (p 3 ~ p^) at the same values 
of x and y. This actual pressure rise may be expressed by the 
coefficient 



/ ^ \ 

^3 *^1 ; = c = actual pressure rise coefficient 

1 .. 2 ^ (for given values of x and y) 

2 P ' V 3 



The difference (c^* ~ c p) is accounted for in part by the gross 
electrical power output as expressed by the previously analyzed 
coefficient c q and in part by additional mixing and friction 

losses that characterize the real system. 

The mass flow rates across sections 1, 2 and 3 may be written 



pA 1 V 1 = xpA 3 V 3 


(4.5) 


pA 2 V 2 = (1-x) A 3 V 3 


(4.6) 


PA 3 V 3 


(4.7) 



From these relations and from Eqs (4.1) and (4.2) we find 
the two area ratios 
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These ratios satisfy the constant area condition 




We can solve these last three equations for the velocity 
ratio 

ll\ = (l-x)y 

V 3y / (y-x) (4.11) 

The momentum equation can now be applied to the ideal ejector 
as follows 

(p 3 * - Pl )A 3 = pA 3 V 3 2 | x (^j * (1 - x>(^) 

2 

We next divide through by ^.PA^V^ and use Eqs 
to eliminate the velocity ratios (V^/V^) and 
simplification, this reduces to the important 



- 1 (4.12) 

(4.2) and (4.11) 

(V 2 /V 3 ) * After 
result 



. * = 2 x(y - l ) 
c p (y - x) 



(4.13) 
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Next consider a purely hypothetical reversible device which 
is not an ejector at all but a type of ideal turbomachine. This 
device receives two streams of fluid identical in all respects to 
the streams at stations 1 and 2 of the ideal ejector. It dis- 
charges a stream identical in all respects to the stream at sta- 
tion 3 of the ideal ejector except for the pressure p^ 3 which 



is higher than p 3 . Again the density p is taken as constant. 
The reversible device operates adiabatically and hence isentropi- 
cally. Under these conditions the following energy equation holds. 
Recall that p^ = p^ • 



The above equation can be simplified by subtracting from it, 
term by term, the following identity. 



Upon carrying out the above subtraction, then dividing through 



find the pressure rise coefficient for the reversible device to 
be 





(4.14) 



xm 3 p 1 + (1 - x)m 3 p x = m 3 p x 



(4.15) 



. 2 

by the reference power m 3 (V 3 /2) and simplifying, we readily 




C. 



PR 




(1 - x) 3 v 2 
(y - x) 2 



(4.16) 
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Comparison of this result with Eq. (4.13) shows that even the 

ideal ejector has a certain inherent mixing loss, call it c * , 

m 

as compared with a truly reversible device. Thus 



c 



* 

m 



(C PR " C p } 



(4.17) 



Of course the mixing loss 

★ 

be somewhat larger than c 

m 

it in the form 




of the actual ejector will 
m J 

and it is advantageous to express 



(4.18) 



where n E / which we term ejector effectiveness, is a dimension- 
less parameter smaller than unity which must be determined from 
experimental data. 

Eq. (4.18) identifies one of the component losses that reduce 

the actual pressure rise coefficient c below its ideal rever- 

P 

sible value c pp . Another loss is that associated with separa- 
tion, friction and secondary flow effects within the ejector; let 
us denote it by symbol c^ E . We have seen that a third component 
is simply the gross electrical power output as expressed by the co- 
efficient c g . Thus the overall energy balance within the ejec- 
tor requires that 



C P C PR (C PR C P } 



c fE " c e 



(4.19) 



The resulting pressure rise through the ejector as expressed 
by Eq. (4.19) must be just great enough to offset the losses in 
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the return passage. These again are losses associated with separa 
tion, friction and secondary flows; let us denote them by coeffi- 
cient c^ R . Accordingly, for the return passage we may write 



c = 
P 



C f R 



(4.20) 



We now eliminate c between Eqs . (4.19) and (4.20) and bring 

P 

all terms to one side. We also introduce the following definition 
for the total friction loss in the fluid flywheel, namely. 



c f c fE + c fR 



(4.21) 



The result is 



'PR 



★ 

— (C pR - C ) + C, 

ru PR P J 



Now using Eqs. (4.13), (4.16) 

components in terms of variables 



c = 0 (4.22) 

e 

and (3.12) to express the above 
x and y , we obtain 



(1 - x) 3 y 2 
(y - x) 2 

. L 

j (1 + 2 X ) 





2 

xy + 



(1 - x 



3 2 

v 



(y 



_ >, 2 
X/ 



1 



2x (y 

“ly 




(4.23) 



This is the overall energy balance equation for the "fluid 
flywheel" and it is fundamental. Notice that the term in the first 
pair of curly brackets represents the ideal power output of a re- 
versible device, the term in the second pair of curly brackets 
represents the sum of the fluid power losses caused by mixing and 
friction and the term in the third pair of curly brackets represents 
the gross electrical power output. 
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conversion 



Eq. (4.23) suggests the utility of defining a 
efficiency" n cv as the ratio of the gross electrical power c g 
to the ideal power c pR . Thus 



n 



cv 




{ v 



(1 + 5x ) 
(1 + 2x) 2 



xy 




(1 - x) 3 y 2 
(y - x) 2 




(4.24) 



Notice that low values of the parameter v , which is ultimately 
determined by the electrical breakdown limit, will necessarily 
produce low conversion efficiencies. 

In connection with Eq. (4.23) and (4.24), notice that x 
always falls between zero and unity while y is always greater 
than unity but has no preassigned upper limit. 

Assuming that parameters ' c f an< 3 v are known or can 

be estimated, Eq . (4.23) then fixes a corresponding unique relation 

between variables x and y . This may be described by a line in 
the xy plane. All subsequent calculations must be confined to 
points that lie along this line. 

If the area ratio _1 be already fixed, then we have from 

A 3 

Eq. (4.8) the additional constraint 




Of course for a fixed value of (Aj/A^) , Eq- (4.25) also 
represents a line in the xy plane. In fact it is a straight 
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line of slope (A^/A^) which passes through the origin. The one 
or more points of intersection of this straight line with the 
curve associated with Eq (4.23) fix the conditions at which steady 
operation of the fluid flywheel is possible. The location of such 
a point determines the values of x and y and hence ultimately 
determines all other performance characteristics of the system as 
well . 

Notice that Eq. (4.23) lends itself to further algebraic sim- 
plification which makes the subsequent numerical solution easier. 
On the other hand the terms in the simplified result cannot be 
given the clear physical interpretation that is possible with 
Eq. (4.23) itself. After simplification, Eq. (4.23) becomes 



1 

Hi 



2x(y - 1) 
(y - x) 



- i - 1 



2 (1 
xy + - — 

(y 



, 3 2 
Y- 



- 1 



- x) 



- V 



(1 + 5x) 
(1 + 2x) 



xy - 



c £ = 0 



(4.261 



For purposes of numerical solution, it is useful to cast 
Eq. (4.26) into a more tractable form as follows. Let 



= 2(y-D 



(4.27) 



Ao = ± - l y 



(4.28) 



a -3 = vy' 



(4.291 
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(4.30) 



A. = c ^ — + 1 

4 f n E 



With this notation, Eq . (4.26) may be rewritten as follows, 

where F is a function that must equal zero. Thus 

3 n 



F = A, , ft , - A- 

i (y-x) 2 



X + 



(1-x) 

(y-x) 



- A ( _a =0 
a 3 ? "4 u 

(l+2x ) 1 4 



(4.31) 



2 2 

Next we multiply through by (1 + 2x) (y-x) and denote 

the resulting function by G. Upon expanding the right side, 
regrouping terms and simplifying, we find that G is a quartic 
in x. The coefficients of this polynomial turn out to be 



B 4 = - 4A x -4 (3-2y) A 2 -5A 3 -4A 4 (4.32) 

B 3 = 4 (y-l)A 1 -4y (y-2)A 2 +(10y-l)A 3 +4 (2y-l)A 4 (4.33) 

B 2 = (4y-l) A^- (4y^-2y-5 ) A 2 ~y (5y-2 ) A 3 ~ ( 4y 2 -8y+l) A 4 (4.34) 

B 1 = Y A 1 -(y 2+1 ) A 2 - ^ 2A 3" 2 y {2 ^ _1)A 4 (4.35) 

B 0 = _A 2' y2A 4 (4,36) 

The function G may then be written 



G = B 4 x 4 + B 3 x 3 +B 2 x“+B^x-!-Bq - 0 



(4.37) 



The unknown root x can be found by Newton's method. For 
this purpose we require also the derivative of G which is 



G' = 4B 4 x 3 +3B 3 x 2 +2B 2 x+B. 



(4.38) 
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n 



Let x be any trial value of x and let G and G' 
n n 

be the corresponding values computed from Eqs. (4.37) and (4.38) 
Then according to Newton's method, the next trial value x n+ j_ 
should be taken as 



G 

„ _ n 

n+1 X n G' 

n 



(4.39) 



This procedure converges rapidly to the true root x . 

This method was used in calculating the numerical example which 
is presented in a later section. 
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5. Thermodynamic Performance 



As shown in the schematic temperature entropy diagram, Fig 
5:1, the tluid flywheel receives a supply m^ (kg/sec) of super- 
heated steam at condition p , T and discharges saturated or 

o o 

nearly saturated liquid at condition p^ , . It gives up heat 

to the ambient atmosphere at absolute temperature T A and delivers 

a gross electrical power output P Q (watts) . It also requires 

an input of electric power, call it P , to energize the injec- 

tor. This is normally a small percentage of the gross electrical 

output and is conveniently expressed by means of an excitation 

efficiency n such that 
x 




(watts ) 



(5.1) 



(continued on page 5-3) 



TEMPERATURE T 




Fig. 5.2. Schematic Temperature Entropy Diagram for Primary 
Stream. 
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After being condensed and leaving the fluid flywheel, the 

liquid water must be pumped back into the boiler. For an ideal 

pump of 100% efficiency this would require the expenditure of shaft 

power in the amount 
* 

P p = m l (p o “ p 4 } v 4 (watts) (5.2) 

where denotes the specific volume of the liquid at state 4. 

For a real pump of efficiency n the corresponding power is 

ir 

* 

p 

p = — 2— 

P Hp (watts) (5.3) 

Thus it is seen that for the overall cycle, the net useful 
power output may be written as 

P e" = (n x P e ' P p* /7 y (watts) (5.4) 

The rate of heat input to the steam passing through the boiler 
and superheater is 

Q = m, (h -h.) - P^ /n (watts) (5.5) 

1 O 4 P P 

where h and h. denote the enthalpies (joule/kg) at states 
O 4 

o and 4 . 

For the hypothetical case of 100% pump efficiency, we denote the 

• ★ 

corresponding heat rate calculated from Eq. (5.5) by symbol Q . 

This differs only very slightly from Q . 
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The overall cycle efficiency now becomes 



n 



c 




- /n T 



(5.6) 



It is of interest to compare the above performance of the 
actual EGD device with the ideal performance of a hypothetical 
reversible device operating under the same input, output and am- 
bient conditions. According to the theory of availability in 
steady flow, the maximum electrical power that could be produced 
by such an ideal device is 



* 



P 

e 






(watts ) 



(5.7) 



where s Q and : s^ denote the entropies at states o and 4, 
respectively (joule/kg °K) . 

The cycle efficiency of this ideal device would be 



★ 




It is also instructive to compare the net output of the real 
system with that of the above ideal and reversible device. We 
term this ratio the relative output. Thus 

rel U*-p* 

\ e p 

All of the foregoing relations can be put into dimensionless 
form by dividing all quantities having the units of power by a 
suitable reference power. In the previous analysis of the fluid 
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flywheel, it was natural and convenient to choose this reference 
as the kinetic power of the total stream at the ejector exit, 
station 3. While it would be possible to retain this same reference 
in this section, it is advantageous to choose here as reference, 
the kinetic power of the primary jet at the ejector inlet, station 
1. This choice reduces the final results to a simpler form and 
makes them somewhat easier to interpret. Thus the reference power 
is now 



m 






h l) 



(5.10) 



In computing this reference power it is permissible to assume 
that state 1 is fixed by the two conditions 



S 1 = S o 



(5.11) 



Pi - P 4 



(5.12) 



In connection with Eq. (5.10) it is also instructive to note 



that 



2 1 2 M _ 
V. "ZT a Mn 2 

,y, H^_l_2ol 

(h o h l 2 



i + 



(5.13) 



Recall also that if we introduce the optimizing condition 
developed earlier, namely, 



2 2 1 
M = M = — 

1 lopt 2 



:5.14) 
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then Eq (5.13) 




simplifies further 




2 (y+3) 



to 



(5.15) 



Incidentally, it is also worth noting that the overall pres- 
sure ratio across the fluid flywheel is 




Moreover, if we observe the restriction of Eq (5.16) this becomes 




For example, for y = 1.3 this gives 



(5.17) 



= 0.731 (5.18) 

This result expresses one of the fundamental limitations of 
the fluid flywheel concept, namely the restriction to a fixed 
pressure ratio. The reason for this restriction was explained 
earlier in connection with Fig. 2.1. 
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Upon simplifying and reducing the various foregoing thermo- 
dynamic relations in the manner explained above, we obtain the 
following results. 



e _ * _ i (T 4 " T a ) ( s o" s 4 } ' Reversible 

= K — { 1 + ■ 



(h -h, ) 
l o 1 



-or^r 



Electric Power 



(5.19) 



m. (h -h ) 
1 o 1 



- = K 



* \ <Pq-P 4 ) ''4 

<V h l> 



Reversible 
Pump Power 



(5.20: 



* * 

Q 



vw 



= K, 






( W K * 

(h -h,) " *P 



Reversible 
Heat Input 



(5.21) 



m, (h -h. ) 
1 o 1 



= K, 



<W fc 

( V h l’ °p 



Actual Heat 
Input 



(5.221 



m (h -h. ) 
1 o 1 



= K = v 



(1+5 x) 
( 1+2 x) 



Gross Electric 
Power 



(5.23) 



m (h -h. ) 
1 o 1 



= K 



s_ 

n xV n 



Net Electric 
Power 



(5.24) 



K 



n c K q \ n x X e 



Cycle Efficiency 25 ^ 



T * * 

-* < K e - K p » 

K Q 



Reversible 
Cycle Efficiency 



(5.26) 



rel 



n k - k /n 
x e P S 
* * 

K — K 
e n p 



Relative Output 
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Once Eq. (4.37) of the previous section is solved, it becomes 
a relatively straight forward matter to apply the relations of 
the present section to calculate the various thermodynamic per- 
formance parameters. Recall that Eq. (4.37) defines a whole series 
of points in the xy plane. By calculating the thermodynamic 
performance parameters of this section for such a series of points, 
and by comparing the respective results, it becomes possible to 
choose that optimum design point which yields the best overall 
trade off among the various desired performance characteristics. 

So far only an initial trial calculation of this type has 
been made on the basis of a plausible initial set of key para- 
meters. This result is summarized in a later section. 
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6 . Summary of Principal Calculation Formulas 



The principal equations developed in the preceding sections 
are summarized below for easy reference. They are listed without 
further explanation in the approximate order in which they would 
be used in a typical calculation. 

Input Parameters 



1 




( 6 . 1 ) 




( 6 . 2 ) 




2 



( 6 . 3 ) 



v = yip 8 



( 6 . 4 ) 



-j- = ( h 0 * h !) (where s 1 = s Q , p x = p 4 ) 



( 6 . 5 ) 



K 



e 




( 6 , 6 ) 



K 



( Pq-P 4 )V 4 



( 6 . 7 ) 



P 




( 6 . 8 ) 



K, 
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k q = 



( V h 4 ) _ ^ 



(h -h. ) n 

o 1 p 



iC* K e - K , 



Solution of Energy Balance Equation 

A. = 

1 n E 



Ao = [ 7T ~ 1] Y 2 



A. = vy' 



A 4 = C f ' ^7 + 1 



B 4 = -4A 1 -4 (3-2y) A 2 -5A 3 -4A 4 



B 3 = 



4 (y-l)A 1 ~4y (y-2) A 2 +(10y-l)A 3 +4 (2y-l)A 4 



B 2 = 



(4y-l) A 1 ~ (4y 2 -2y-5) A 2 ~y (5y-2) A 3 ~ (4y 2 -8y+l) A 4 



B 1 = 



yA 1 -(y 2 +l)A 2 -y 2 A 3 -2y (2v-l)A 4 



B o = 



-A 9 -y A 



G = 
G' = 



B^^+B-^x^+B-pc 2 +B^x+B q 
4B 4 x 3 ^3B 3X 2 +2B 2 x+B 1 



x n+l 



'n 



'n 



n 



(Iterate until x converges) 



(6.9) 

( 6 . 10 ) 



( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 

(6.14) 

(6.15) 

(6.16) 

(6.17) 

(6.18) 

(6.19) 

( 6 . 20 ) 
( 6 . 21 ) 
( 6 . 22 ) 
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Output Parameters 



* 



II 

K 

e 


= < 


n x v “♦**> 
(1 +2x) 2 


K 

_ P 

n p 


Net Electric 
Power 


(6.23) 


n c 


= 


K 

e 

k q 




Cycle Efficiency 


(6.24) 






K " 
e 




Relative Output 


(6.25) 


n rel 




★ * 
< K e - K p » 




K 

e 


= \ 


|v (1+5X) ,| 
! (1+2x) 2 J 




Gross Electric 
Power 


(6.26) 


k pr 


= ' 


+ (1-x) 3 
x ty-x) 


1 | 

2 , 
xy 1 


Ideal Pressure 
) Rise 


(6.27) 


^cv 


- 


K 

e 

k pr 




Conversion 

Efficiency 


(6.28) 



Note: Power coefficients denoted by symbol K are expressed 

in terms of the kinetic power at station 1. Those denoted by 
symbol c are expressed in terms of the kinetic power at station 
3. The conversion between these two forms is simply 

K = -^2 (6.29) 

xy 



6-3 



7. Some Typical Experimental Data and Preliminary Numerical 



Results 



Reference (1) reports measurements of breakdown field 
strength as a function of pressure at constant temperature and 
electrode separation. Their results show that = 9 x 10^ v/m 

and Cg = 9.5 x 10 2 m 2 °K/C . Note the interesting fact that 
steam and air appear to have identical breakdown properties ! 

A system now under development by Marks Polarized Corporation 
is intended to operate with steam at about 100 atmospheres pres- 
sure and sufficient superheat to give about 5% moisture at state 
1. The data in Table 7.1 suggest that under these conditions, 
approximately, from Eqs . (6.1) and (6.2), 



1 




0.920 x 10°_ 




591 




1.273] 


0.953 x 10 4 




10. 


34 x 10* 



0.0070 ~ 0 



(7.1) 



8 




fa. 85 x io -12 ] fo. 


953 


*_10_ 4 ] 




LiO. 34 x 10 6 


1 


[59 i] 


I 2 



= 0.0238 



(7.2) 
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Then from Eq. (6.3) we obtain 






3 

Y 





( 1 + 0) 2 



= 1.32 



(7.3) 



In the absence of adequate test date, the best estimate of 
y that is possible at this time is simply 

y as 1 (7.4) 

According to Eq. (6.4) this gives the estimate 

v = yi^B = (1) (1.32) (0.0238) = 0.032 (7.5) 

The value of the flow loss factor c^ for the annulus may 
be estimated from the data in Ref. (2) on losses in pipe bends. 

At a Reynolds number of 200,000 or above, these losses are mainly 
a function of R/r, where R is bend radius and r is pipe 
radius. For the torus in Marks' design R/r = 5.7. Ref. (2) 
shows that at this value of R/r the loss factor for a 180° 
pipe bend is about 0.4. Thus for a 360° turn the loss factor 
would presumably be twice the value or approximately 

c f » 0.8 (7.6) 

There are qualitative grounds for hoping that the loss in 
a true torus might be somewhat smaller than the above pipe bend 
data would suggest, but there is no hard data available on this 
point. Thus Eq. (7.6) represents a plausible and conservative 
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working hypothesis that can be advanced at this time. 

Values of rig can be estimated from available test data 
on ejector performance. See, for example. Ref. (3). These 
data have not yet been reviewed critically. All we can say at 
present on the basis of general engineering judgment and experi- 
ence is that for the range of ejector design parameters of in- 
terest in the present context, is unlikely to lie above 

0.9 or below 0.5. It is thought reasonable as a tentative ini- 
tial working hypothesis to choose 

n E a 0.8 (7.7) 

Also required for the calculation are certain thermodynamic 
properties. These are summarized in Table 7.1. These data hap- 
pen to be in English units but this is permissible because the 
final results, being dimensionless, are independent of the system 
of units used. For some purposes the steam was treated approxi- 
mately as a perfect gas with a ratio of specific heats v = 1.3. 

The information in Table 7.1, along with the formulas of 
Section 6, permit the various key input parameters to be esti- 
mated as summarized in Table 7.2. These values are regarded as 
the best estimates that are possible at the present time. These 
input data suffice to permit a corresponding performance calcula- 
tion to be carried out. This has been done and the main results 
are summarized in Table 7.3 and in Fig. 7.1. Notice that each 
row or Table 7.3 corresponds to a different value of the ejector 
area ratio A^/A^. 
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Table 7.1 Thermodynamic Properties 



State 




0 


1 


4 


Ambient 


P 


psia 


1500 


1100 


1100 


— 


T 


°R 


1063.9 


1016 


1016 


530 


h 


Btu/ lbm 


1180 . 7 


1156.3 


557.4 


— 


s 


Btu/lbm°R 


1.3470 


1.3470 


0.7575 


— 


P 


lbm/ft 3 


— 


— 


45.45 


— 



Table 7.2 Estimated Values of Key Input Parameters 



Y 


— 


1.3 


n p 

★ 


— 


0.9 


y 




1 


K 

e 

* 




12.7 


V 


~ 


0.032 


K 

P 

* 




0.067 


n E 


— 


0 . 8 


k q 


= 


25.5 


C f 


= 


0.8 


k q 

* 


= 


25 . 5 




= 


1 




= 


0.50 
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Table 7.3 Estimated Performance of EHD Generator 
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The curve ABCDEF shown in Fig. 7.1 is the locus of points 
x,y which satisfy the basic energy balance relation of the 
ejector and fluid flywheel, that is, Eq. (4.23). In view of 
the tentative nature of the values assigned to parameters v, 
n £ and c f , the accuracy of this curve is somewhat uncertain, 
but its general qualitative features are undoubtedly correct. 

Any straight line through the origin, such as line OEB , for 
example, represents Eq. (4.25), and the inverse of its slope 
represents the area ratio A 1 / A 3 • Notice that such a straight 
line may in general intersect the curve at two points such as 
E and B , for example. It can be shown that only the lower 
point E represents stable equilibrium; point B represents 
unstable equilibrium and should be disregarded. For this reason 
the entire dashed portion ABC of the curve must also be dis- 
carded. Only the solid portion CDEF has any practical signifi- 
cance . 

Notice that the tangent point C represents the lowest value 

of A^/A^ and the highest value of y that can be achieved in 

steady state operation. Recall that the gross electric power is 

2 

proportional to the factor ( l+5x) / ( l+2x ) . It can be shown that 

this factor reaches its peak value of 1.04 at x = 0.1 . It is 
therefore desirable to choose the operating point on the curve 
as close to x = 0.1 as possible. In this case the nearest point 
is at D where x = 0*36. Note that this is the lowest value of 
x that can be attained for the assumed values of v , n E and c^ . 
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The corresponding value of the factor ( 1+5 x) / ( l+2x) ^ is 0.95. 

Notice from the data in Table 7.2 that an ideal reversible 
device operating on the given basic thermodynamic cycle would 
produce a dimensionless net power output of amount 

(K e * ~ Kp*) = 12.7 - 0.07 as 12.6 (7.8) 

and would produce an ideal cycle efficiency of amount 

n c * = 0.50 (7.9) 

These theoretical performance limits of the ideal reversible 
device are in themselves very favorable. Unfortunately, the per- 
formance of the actual EHD system as summarized in Table 7.3 is 
seen to be very unsatisfactory. In fact the net power outputs 
and cycle efficiencies are negative for all possible operating 
points! This means that even at 100 atmospheres pressure, the 
phenomenon of electrical breakdown so severely limits the gross 
electrical power output that it it insufficient even to meet the 
demand for power to drive the pump, let alone provide any useful 
output . 

It might be supposed that a radical increase in the electrical 
breakdown strength, if it could be achieved, would in itself trans- 
form this pessimistic picture into an optimistic one, but this is 
not entirely the case, at least not for the present EHD scheme. 
Suppose for example, that the breakdown parameter v could be in- 
creased by a factor of roughly 30 so as to bring the gross elec- 
trical power coefficient up say to 

K e * 1 (7. 10) 
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In that case the cycle efficiency would still be only about 




(7. ID 



or less than 4% . The main reason for this unsatisfactory result 
is that the present cycle squanders available energy by dumping 
heat from the condenser at high temperature to the ambient air 
at low temperature without extracting any useful power from it. 

The irreversibilities inherent in the ejector also waste available 
energy. Moreover, friction losses in the fluid flywheel appear 
to be disappointingly high. As a reflection of these losses, 
note from Table 7.3 that the conversion efficiency n cv remains 
below 4% for all possible stable operating points. 

It is also worthwhile at this point to review the various 
factors that affect the important parameter v which so severely 
limits the electrical output. These factors are revealed clearly 
in Eqs . (6.2), (6.3) and (6.4). It is convenient to combine them 

here in the following wav, assuming ct - 0. Thus 



As y varies from 1 to 1.67, the function of y i n the 
curly brackets above varies from 1.82 to 0.97. Hence there is 
no possibility of radically changing the order of magnitude of 
v by manipulating y . 




(7. 12) 
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Similar considerations apply to the quantity e . It has 
been suggested, see Re. 4, that e might be increased by using 
a polarizable aerosol. But the overall increase that can be 
achieved in this way is expected to be less than a factor of ten 
and, as we have seen, this in itself is not enough to overcome 
fully the limitations of the present EHD scheme. 

An obvious way to increase v is to increase p Q . We have 
noted, however, that even at 100 atmospheres pressure, v is 
still much too small. Any increase beyond this value is not 
feasible for two reasons. Firstly, structural and practical 
problems become troublesome at such high pressures. Secondly 
the linear relation between electrical breakdown strength and 
fluid density as given by Eq. (2.5) breaks down near pressures 
of 100 atmospheres. See Ref. 1. 

The two remaining factors at our disposal are Cg and T Q . 

2 

Note that v is proportional to the quadractic factor (C_/T_) . 

The breakdown constant C is a property of the medium. The 
question whether any aerosol mixture exists that has a substan- 
tially higher value of Cg has not been investigated. It is 
recommended that such an investigation be carried out. 

Finally consider the factor T Q . In the present scheme, 
charged particles of low mobility are created by condensation. 
This requires that T Q be high enough to produce about 5% 
liquid at state 1. This restriction on T Q could be eliminated 
by abandoning the condensation method and utilizing instead the 
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direct injection of liquid droplets or solid particles. This 

change also entails changing the carrier fluid from steam to 

air or some other noncondensing gas . Note that a decrease in 

T from 1016 °R to 530 °R would increase v by a factor of 
o 

about 3.7. 

It seems possible that some combination of the above methods 
could perhaps increase v by a factor of say 20 or so. We have 
seen that this m itself would not make the present fluid fly- 
wheel configuration successful. Nevertheless, an improvement of 
this order might very well suffice to make some alternative EHD 
scheme practical. 
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8. Conclusions and Recommendations. 



Broadly speaking, one of the inherent limitations of EHD 
power generation is the severe restriction imposed by the 
phenomenon of electrical breakdown. Because of this effect, 
the electrical work that can be extracted per unit mass of work- 
ing fluid is very small. In general there are two distinct ap- 
proaches for attempting to overcome this limitation. One method 
is based on multi-staging, that is, on passing the same working 
fluid through a whole series of successive electrical conversion 
sections. The other method is to adopt an ejector system which 
utilizes a primary stream of high velocity but low mass flow 
rate to drive a secondary stream of low velocity but high mass 
flow rate. The conversion of flow energy to electric work is 
accomplished in the resulting low velocity stream. Although 
electrical work per unit mass of low velocity fluid remains low, 
the electrical work per unit mass of high velocity primary 
fluid may be boosted to a sufficiently high value. The effective- 
ness of this scheme is limited to some extent by the character- 
istic ejector mixing losses. 

In this report, a proposed type of single stage EHD genera- 
tor is analyzed which utilized the ejector principle. However, 
this design combines the ejector and the electrical conversion 
section into a single unit; in particular, charged particles are 
introduced into the primary stream only and are subsequently 
collected at the ejector exit section. 

The analysis shows that, unfortunately, this type of design 
gives extremely poor performance. In fact che electrical power 
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generated is insufficient to drive the feedwater pump, let alone 
produce any net useful output. This poor result can be attribu- 
ted to the fact that the charged particles are introduced into 
the primary stream, thus generating a very intense local 
electrical field at this point. 

An obvious remedy for this difficulty is to place the elec- 
trical conversion section after the ejector where the electrical 
field can be much lower. The Appendix evaluates this revised 
design. It is shown in the Appendix that under the most opti- 
mistic assumptions, the cycle efficiency is less than 4% which 
may be insufficient to compete with alternative power systems. 

It may be possible to multi-stage the fluid flywheel 
generators and thereby show an order of magnitude improvement 
in performance. It is recommended that systematic parametric 
studies be carried out to determine the performance possi- 
bilities of such a multi-stage system and to estimate optimum 
design values of the key parameters. It would also be ad- 
visable in the interest of orderly development to restrict this 
proposed next step to the same basic medium, namely steam, and 
to the same basic Rankine cycle as used in the previous studies. 
Of course there is no assurance that a multi-stage design will 
necessarily prove successful because, even if the theoretical 
cycle performance is favorable, practical considerations of 
bulk and complexity might nullify the theoretical advantages. 
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APPENDIX: Performance Characteristics of Revised EHD Generator 



The analysis in the main text clearly indicates that intro- 
ducing charges in the primary stream only, limits the thermody- 
namic state 1 to a relatively high fixed pressure ratio 
(P]_/P Q = 0*731 for = 0.71). An obvious improvement is to 
uncouple the ejector from the electrical conversion section. 

This is done in the following analysis in which the electrical 
conversion section is located after the ejector. 

The work in this Appendix is based on a fixed upper pressure 
p Q = 1500 psia and a fixed steam quality of 0.950 at station 1 
(See Fig. 5.1). However, the back pressure p^ is allowed to take 
on a large number of possible values ranging from 1400 psia 
down to 14.7 psia. At any given back pressure, the area ratio 
A^/A^ is allowed to vary as required to obtain best performance. 
It is found that best performance corresponds to the minimum 
possible value of the mass flow ratio (m^/m^) = x. This is 
governed by the energy balance equation for the fluid flywheel. 

Three cases are considered, as listed in more detail on 
page 2A. Case A represents the present best estimate of the 
three key parameters n^ R , and Cg. Case B represents a much 
more optimistic estimate of parameters n R and c^, perhaps values 
which can never be achieved practically. Case C retains these 
optimistic values of n R and Cf and assumes in addition that the 
electrical breakdown constant Cg can be increased by a factor 
of /TO. Thus, case C represents optimism plus a fundamental 
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REVISED EHD ANALYSIS 



BASIS 



* Design change: electrical power section now placed after 

ejector. 

* Rankine cycle with - 

P q = 1500 psia 

steam quality at state 1 = 0.95 
OPTIMIZING PARAMETERS 

* Pressure p^ in fluid flywheel 

* Ejector area ratio A^/A^ 

CASES STUDIED 







A 


B C 


Ejector Effectiveness 


n E 


0.8 


0.9 0.9 


Friction Factor 


C f 


0.8 


r— 1 
O 

» — 1 
O 


Electrical Breakdown 
_ Constant 
(irT °K/Coulomb) 


C B 


9. 49xl0 3 9. 


49xl0 3 /l0x9.49xl0 3 



2A 



scientific "break-through". 

The single most important result of the analysis is the 

curve of cycle efficiency r) versus back pressure p. , Fig. 1. 

c _L 

This result confirms that the proposed design revision has indeed 
improved the performance over that of the original version. Re- 
call that the original design yielded a negative cycle efficiency! 
The revised design, case A, produces positive efficiencies for 
values of p^ above about 900 psia. The peak efficiency occurs 
approximately at p., = 1400 psia. Unfortunately, this peak value 
is still much too low, below 0.003. As expected, case B shows a 
marked improvement over case A, but nevertheless n c still re- 
mains below 0.01. Even the extremely optimistic case C gives a 
peak cycle efficiency of less than 0.04. 

The efficiency curve for case C has an interesting bimodal 
form with the solution above p^ = 1000 psia quite different from 
that below this limit. The analysis shows that the ejector cannot 
function for case C for values of p^ above about 1000 psia, so 
the solution over this region pertains to a simple EHD conversion 
in a single stream; the ejector principle cannot be utilized in 
this region. It is ironic that the highest efficiency shown, 
nearly 0.04, occurs at about p^ ~ 1400 psia where the fluid fly- 
wheel concept does not even apply! A second maximum cycle effi- 
ciency of about 0.026 occurs at p^ = 270 psia. This point repre- 
sents normal operation which utilizes the fluid flywheel principle 
to good advantage. 
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FIG IA EFFICIENCY VERSUS PRESSURE AT STATION I. 
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The results show that the final cycle efficiency that can 
be acheived, and the optimum back pressure p^ and ejector area 
ratio (A^/A^ ) that will produce this optimum depend critically 
upon the three key parameters n E , c^ and Cg. If these para- 
meters can be estimated accurately then design optimization and 
resulting performance can be accurately predicted and controlled, 
otherwise not. 

The results also show, however, that even under the most 
optimistic assumptions regarding n £ , c^ and C , the cycle 
efficiency is disappointingly low - less than 0.04! It is 
therefore clear that, viewed as a single stage device, the pro- 
posed EHD generator, even with the noted design improvements, 
has performance limitations which are just too severe to permit 
it to compete successfully with alternative power systems. 

On the other hand, the possibility still remains that a 
multi-stage version of this generator might show a further order 
of magnitude improvement in perf ormance . It is not clear at 
this time whether the ejector and fluid flywheel concepts can 
be incorporated to advantage in such a multi-stage design. 

The remainder of this Appendix summarizes the analysis of 
the revised scheme. The notation and approach is generally 
similar to the main text. The derivations, formulas, tabula- 
tions and diagrams are complete but no attempt has been made 
to integrate them into a smooth flowing text. 
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ELECTRIC POWER COEFFICIENTS 



P e 2 V 1 

5^- = c a = vy where y= (^=-) 



■K3 



pT- = K e = ^ where x = (— ) 
KI 



v = 0.5837- 



( h o - h.: 



mi , 
m3 

(CASES A & 3) 



= density, lbm/ft^ 

(h Q - h^) = enthalpy drop, Btu/lbm 



Note: Above constant 0.5837 is based on a value of the break- 

down constant of: 

Cg = 9.49 x 10^ m^°K/coulomb 

An assumed increase in C g by a factor v'To would increase v by 
a factor of 10. This would give: 



v 



5.837- 



(h 



- h x ) 



(CASE C) 
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ENERGY BALANCE 



OPERATING LIMIT 



The fluid flywheel can function only if parameters n E , c^ 
and v satisfy the inequality 



c-(-i - 1) + v (-i + 1 + c-)<l 
f Rp n- f 



(l) 



OPERATING RANGE 



Let 



C 7 = (-=■ - 1 + v) 

2 n E 



C 1 = -i 

n E 



( 2 ) 

(3) 



C = (-± + 1 + c-) 



(4) 



2 2 

R = C7 - 4C-,C 
1 2 o 



(5) 



A solution exists for any valve of y between the limits - 



C 1 - R 
■^min ^ 2C 2 ^ 



C. + R 

y = ( 

^max 2C 2 



( 6 ) 

(7) 



MASS RATIO x! 

The solution for mass ratio x is found by the following 
sequence of calculations. 



, 2 (y. 

l 



A 2 = (-p - l)y 2 
2 n E 



( 8 ) 

(9) 



A-, = (- — + 1 + Cr + uy 2 ) 



( 10 ) 
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( 11 ) 



B 2 = A x + (3 - 2y)A 2 + A 3 



B 1 


- ^ A 1 ■ 


- (y 2 - 3)A 2 + 2yA 3 


(12) 


B 

o 


= A 2 t 


y 2 A 3 


(13) 


s 2 


= *?- 


4B,B 
2 o 


(14) 




B, ' 


- S 




X 


= ( ^ 
v 2B. 


) 

3 


(15) 



Best performance corresponds to the minimum possible value 
of x. This is best found by numerical trial and error. 
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THERMODYNAMIC PARAMETERS 



Let 




= power, watts 

= kinetic power at station 1, watts 



Then 



K 



K2 



= power coefficient 



K = {1 + 
e 



<T 4 - V (s o - s < 

(h o - V 



* <Pp - ?4 )v . 

p (h - h, ) 



Reversible Electric Power 



Reversible Pump Power 



(h o " h 4> 
K = f Q 

k q 1 (h - h, ) 



K 



Heat Input 



P" 


n v k* 




e _ 


{ ir - ^ }<h o - h i ! 


Specific Output Watt Sec/Kg 


m l 


P 

n ,'■> k* . 




n = 


{ 6 - - —Hi- i. 


Cycle Efficiency 


c 


x V k q 



PRINCIPAL PERFORMANCE PARAMETERS 



Values used in calculations- 

n x v 

Excitation Factor — — = 1.00 



Pump Efficiency 



0.90 



THERMODYNAMIC PROPERTIES 
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THERMODYNAMIC COEFFICIENTS 



Pi 


V 


K* 


psia 


1400 


0.4063 


54.94 


1300 


0.1721 


26 .53 


1200 


0.1001 


17.65 


1100 


0.06274 


12.75 


1000 


0 . 04244 


10.04 


900 


0.02958 


8.151 


800 


0.02073 


6.769 


700 


0.01455 


5.716 


600 


0.01010 


4.876 


500 


0.006839 


4.178 


400 


0.004430 


3.574 


300 


0.002663 


3.037 


200 


0.001391 


2.537 


100 


0.0005091 


2.016 


14.696 


0 . 00004396 


1.374 


Notes : 


Ah 


- (h o - V 




AH 


■ <h o - Vp, 


= 14.67 



k q 


K* 

P 


V AH ; 


110.1 


0.0855 


0.009590 


53.09 


0.0778 


0.2071 


35.30 


0.0732 


0.03241 


25.47 


0.0667 


0.04680 


19.97 


0.0617 


0.06214 


16.28 


0.0568 


0.07940 


13.53 


0.0521 


0.09973 


11.43 


0.0472 


0.1233 


9.762 


0.0424 


0.1513 


8.375 


0.0377 


0.1855 


7.179 


0.0329 


0. 2288 


6.113 


0.0281 


0.2871 


5.124 


0.0229 


0. 3702 


4.083 


0.0169 


0.5211 


2.726 


0.0088 


1.0000 



Since system is closed, atmospheric pressure p 1 = 14.696 
psia has no particular significance but serves as a ^convenient , 
if arbitrary lower pressure limit for the cycle calculations. 
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RESULTS OF EJECTOR OPTIMIZATION 
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PRINCIPAL RESULTS 



pi 

psia 


CASE A 


CASE B 


CASE 


c 


n c 


KW SEC 
Kg 


n c 


P"/ih. 
e' 1 

KW SEC 
Kg 


n c 


p ;/*i 


1400 


0.0028 


3.61 


0.0030 


3.83 


0.0360 


46.14 


1300 


0.0022 


2.86 


0.0041 


5.42 


0.0308 


41.05 


1200 


0.0018 


2.45 


0 . 0051 


7.06 


0.0260 


36 . 14 


1100 


0.0013 


1.88 


0.0062 


9.05 


0.0217 


31.40 


1000 


0 . 0007 


1.09 


0.0072 


10.85 


0.0183 


27.59 


900 


0.0001 


0.15 


0.0080 


12.61 


0.0186 


29.21 


800 


-0.0006 


-0.96 


0.0087 


14.26 


0.0194 


31.82 


700 






0.0092 


15.70 


0.0206 


35.22 


600 






0.0094 


16.81 


0.0220 


39.35 


500 






0.0093 


17.60 


0.0234 


44.17 


400 






0.0086 


- 17.17 


0.0248 


49.50 


300 






0.0074 


15.69 


0.0257 


54.79 


200 






0.0053 


12.14 


0.0253 


58.27 


100 






0.0018 


4.73 


0.0211 


54.57 


14.696 






-0.0067 


-8.15 


0.0050 


16.65 
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Fig. 2A v versus pressure at station 1 



14A 



1400 




15 A 



Fig. 3 A Enthalpy Drop Relative to that at p = 14.696 





Fig. 4A Typical y versus x 
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Pi 

psia 





opt 



Fig. 5A Area Ratio versus 
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Pi 

psia 




Fig. 6A Velocity Ratio versus 
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Pi 

PSJG 




Fig. 7A Minimum Mass Ratio versus 
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Fig. 8A Net Electric Power versus 
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ELECTRIC POWER 



1 2 2 2 2 2 

P e feC/R^^ = £ C B R^P 3 

P K3 2 P 3 V 3 A 3 V 3 



Also 



V ] 

p 3 = Pl and V 3 = - 



2 „ 6C B 2R2 fl eC B 2r2 <=1 

. ■ c e - vy where v - — j 3 ( h " — h ] 

v . o 1 



Reqrouping factors - 
2 2 

£C r R ( 
V . ( R ) 



(ho V 



(dimensionless ) 



In metric units - 
2 2 

£C R R 1 

; — 2 ) = ±[8.854 x 10 



vm 



-] [ 9 . 4 9 x 10 



3 m 2 °K 
cmb 



x 4 61,- 



xg 



= 84.73 ( 



Jm . 

k 2 ' 
g 



In English units - 

sC 2 r 2 j 3 

( — |— - ) = [84. 73j^j][9. 480 x 10 " 4 

x [3.281— ] 3 [0.4536 T 2 -] 2 
L m J L lbm J 



= 0.5837 



Btu ft~ 
lbm 



Hence 



v = 0.5837 



< h o - h x ) 



(dimensionless) 



where 



where 


3 

= density, lbm/ft 


(h o ' 


- h^) = enthalpy drop, Btu/lbm 


also 


P c 

_e = K = — = 

P K1 6 xy 2 x 


SUMMARY 




V = 0. 


P, p, lbm/ft 2 

c o D 7 -L f -L 1 

(h - h, ) 1 (h - h, ) Btu/lbm ' 

O I o 1 

p ~ 

e c ,,2 

p = e = vy 

K3 

P e _ K - u 
p ki e y 
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ENERGY BALANCE 



By a simple revision of previous work we obtain the two key 
equations . 



(— ) = (-) 
A3 ' V 



( 1 ) 






} 



- {c f + vy z } = 0 



( 2 ) 



Eq (2) can be reduced to the form: 

u. Lk ~ 3 

i 'y - x' 2 



A, { 



X } . A { x + (1 A, = 0 



(y - x) 



(3) 



where 



A = 2 (y - 1) 
1 



a 2 - <4 - i)y2 

A, = (- — - + 1 + C£ + vy 2 ) 



(4) 

(5) 

( 6 ) 



Eq (3) can be further reduced to the form: 



B 0 x - B,x + B = 0 
2 1 o 



3 2 = A^ + (3 - 2y)A 2 + A 3 



(7) 

(3) 



3^ = yA 1 - (y - 3)A 2 



2yA- 



(9) 



s n = A o + y A - 



( 10 ) 



Eq (7) may be solved by the quadratic formula. The larger 
of the two roots is extraneous and may be dropped. 



Let S 2 = B. 2 - 4B-B 
1 2 o 

B 1 - S 

Then X = (— ) 

2B_ 



( 11 ) 



( 12 ) 
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For the limiting case x = 1 , Eq (2) reduces to the form: 



where 



Let 



c 2 y - c lY + c o - 0 



C, = (-± - 1 + v) 



C, = 



C o * + 1 + 



2 2 

R = Cl - 4C~C 
1 2 o 



(13) 

(14) 

(15) 

(16) 
(17) 



The corresponding roots are: 



C, - R 

V min = (— kr ) 

mm C 

C, + R 



y max ^ 2C. 



-) 



(18) 

(19) 



The above roots are real only if R z is not negative. 
This will be the case if and only if : 



c,(- - 1) + v(— + 1 + c-) <1 (20) 

f n E n E f - 
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